Relation between shear parameter and Reynolds number in statistically stationary 

turbulent shear flows 
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Studies of the relation between the shear parameter S* and the Reynolds number Re are presented 
for a nearly homogeneous and statistically stationary turbulent shear flow. The parametric investi- 
gations are in line with a generalized perspective on the return to local isotropy in shear flows that 
was outlined recently [Schumacher, Sreenivasan and Yeung, Phys. Fluids 15, 84 (2003)]. Therefore, 
two parameters, the constant shear rate S and the level of initial turbulent fluctuations as prescribed 
by an energy injection rate €i n , are varied systematically. The investigations suggest that the shear 
parameter levels off for larger Reynolds numbers which is supported by dimensional arguments. It 
is found that the skewness of the transverse derivative shows a different decay behavior with respect 
to Reynolds number when the sequence of simulation runs follows different pathways across the 
two-parameter plane. The study can shed new light on different interpretations of the decay of odd 
order moments in high-Reynolds number experiments. 

PACS numbers: 47.27.Ak, 47.27.Jv, 47.27.Nz 



INTRODUCTION 

Homogeneous shear turbulence can be considered as 
the first non-trivial extension of homogeneous, locally 
isotropic turbulence. It is a flow that is thought as 
a bridge between the strongly idealized homogeneous 
isotropic turbulence and more realistic turbulent shear 
flows such as channel flows [1]. A constant mean 
shear rate S is present which is a large-scale source of 
anisotropy and the question that immediately arises is 
how the statistical properties at the smallest scales of 
the turbulent flow are affected by its presence. After 
the pioneering work of Lumley [2] in which he predicted 
by dimensional arguments a rather rapid decay of such 
anisotropy, namely with ~ R^ 1 , this particular flow came 
under renewed interest in the last decade. Several sys- 
tematic measurements in simple shear flows [3, 4, 5, 6] 
and in atmospheric boundary layers for the largest acces- 
sible Taylor microscale Reynolds numbers of R\ ~ 10 4 
were presented. [7] All experiments detected systematic 
deviations from Kolmogorov's concept of local isotropy. 
[8] The decay of odd order normalized transverse deriva- 
tive moments with respect to Taylor microscale Reynolds 
number occurs with a larger exponent than -1 and the 
co-spectrum of the shear stress deviates from the classi- 
cal -7/3 law. [7] For moderate Reynolds numbers, such a 
persistence in the decay of odd order moments was also 
found within direct numerical simulations (DNS) and its 
connection to coherent structures and intcrmittency cor- 
rections could be addressed. [9, 10, 11, 12, 13] 

Recently, the record of experimental and DNS data on 
homogeneous and nearly homogeneous shear flows was 
collected and discussed anew from a generalized perspec- 
tive by taking into account the role of small-scale in- 
tcrmittency and mean shear. [14] It was found that the 
operating points of all those experiments are scattered in 



a two-parameter plane which is spanned by the Reynolds 
number Re and the shear parameter S* . Different experi- 
ments followed different pathways across such plane caus- 
ing, e.g., variations in the decay behavior of odd order 
moments when simply projected onto the Reynolds num- 
ber axis as it is done usually when the issue of isotropy 
is discussed. This was identified as one possible reason 
for different interpretations of the data in terms of the 
return to local isotropy for higher Reynolds numbers. 

One outcome of Ref. 14 is the necessity for more sys- 
tematic numerical experiments which will be presented 
in the following. Here, two system parameters that de- 
termine the homogeneous shear flow will be varied: the 
constant shear rate, S 1 , and an energy injection rate, 6j n , 
that prescribes the amount of initially isotropic turbulent 
fluctuations. The latter parameter can also be thought 
of as a substitute for active or passive grids that are of- 
ten placed in wind tunnel experiments before the working 
fluid enters the shear straightener (see e.g. Refs. 4 to 6 
and the sketch in Fig. 1). Such additional device (which 
is optional, but frequently used) increases the turbulent 
fluctuations and thus the Reynolds number while oper- 
ating with the same mean shear rate S. Beside this main 
motivation coming from the experiments in nearly homo- 
geneous shear flows, it allows here for reaching a stable 
turbulent regime for small shear rates. Consequently, the 
numerical experiments will give us hints on the following 
functional dependencies 



Re = Fi(u,S,e in ) , 
S* = F 2 (v,S,e in ) , 



(1) 

(2) 



and would therefore allow for a more systematic study 
of the deviations from local isotropy, e.g. how trans- 
verse derivative moments decay along particular path- 
ways across the two-parameter plane. We can simplify 
the functional dependencies in relations (1) and (2) by 
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FIG. 1: Sketch of a typical experimental situation in a sim- 
ple shear flow and its relation to the present investigation. A 
shear straightener sets up a mean flow with a shear rate S and 
an active or passive grid produces additional turbulent fluctu- 
ations. The turbulence is decaying along the downstream di- 
rection. The small volume would be the one we model within 
the DNS. 



keeping the kinematic viscosity v constant throughout 
the study. 

Numerical models for homogeneous shear turbulence 
can be divided into three different groups. The first one 
uses finite difference methods with shear periodic bound- 
ary conditions on the grid in physical space. [15] The other 
two approaches use pseudospectral techniques with pe- 
riodic boundary conditions. Rogallo suggested a time- 
dependent coordinate transformation for the inclusion of 
constant shear in such a simulation of turbulence. [16, 17] 
The consequence is that the numerical grid gets steadily 
skewed and has to be remeshed after time t = 2/S. In 
this method, it was hard to reach a stationary turbu- 
lent state, the integral length scale grows exponentially 
in time until the order of the box length is reached. Sta- 
tistical investigations that continued afterwards had to 
be run for very long time intervals due to large fluctua- 
tions of energy or enstrophy. [9, 12, 18] In Ref. 9 inte- 
grations up to St ~ 10 3 were performed. For very large 
S aliasing errors were reported. [19] A method to over- 
come some of the problems was suggested recently. [10] 
The method avoids the rcmeshing procedure. A statisti- 
cally stationary state with smaller fluctuations of energy 
and enstrophy can be attained by an appropriate vol- 
ume forcing in combination with free-slip boundary con- 
ditions. The boundary conditions will cause small devi- 
ations from transverse homogeneity, so strictly speaking 
our system is a nearly homogeneous shear flow, as it is 
the case for the measurements. 

The outline of the manuscript is as follows. In the next 
section, we present in brief the numerical model and the 
volume forcing schemes which are related to both outer 
parameters, S and e m . In section III we conduct de- 
tailed studies of the statistical stationarity and point to 
differences in comparison with the classical remeshing ap- 
proach. Afterwards, we are going to discuss the relation 
between shear parameter and Reynolds number. Predic- 



tions that are made on the basis of dimensional argu- 
ments will be compared with numerical findings. Based 
on these results, the behavior of derivative moments with 
respect to Reynolds number and shear parameter will be 
studied. Finally, a summary and an outlook are given. 



NUMERICAL MODEL 
Equations and boundary conditions 

With length scales measured in units of the box width 
L y , and time scales in units of the eddy turnover time 
L y /v rmsi the dimensionless form of the equations for an 
incompressible Navier-Stokes fluid become 



^ + (u.V)u = -V P+ ^V 2 u + F 



Vu 



0, 



(3) 
(4) 



where p(x, t) is the pressure, u(x, t) the velocity field. 
The large scale Reynolds number is then 



Re = 



V rmS Ly 



(5) 



The velocity components are decomposed in a mean frac- 
tion and a fluctuating turbulent part, Ui = (u,) + Vi for 
i = x,y and z, the so-called Reynolds decomposition. 
The mean flow profiles are given for the homogeneous 
shear flow by 



(u x ) 



A,T 



Sy, (u y ) 



A,T 



(u z )a,t = 0, 



(6) 



where x are streamwise (or downstream), y shear (or 
wall-normal), and z spanwise directions, respectively. 
Statistical averages are denoted by parentheses for the 
following and indices indicate whether a time average, 
(•)t, a volume average, (-)v, an average over x — z 
planes at fixed y, (-)a, or combinations of them are 
taken. Consequently, the root mean square velocity reads 
Vrms = \J (v 2 )y,T and the dimensionless shear parameter 
becomes 



e„2 

_ v rms 



(7) 



with the energy dissipation rate e = (e(x,t))vT and 
e(x,t) = W2) £ij(a,-«i(x,t) +0 iWj (x,t)) 2 . 

The pseudospectral method is applied. The equa- 
tions are integrated by a second order predictor-corrector 
scheme where the time stepping satisfies the Courant- 
Fricdrichs-Levy criterion. The Courant number was al- 
ways below 1/2. The dissipative term was included as an 
integrating factor exp(— vk 2 t) for every mode Uk(t). De- 
aliasing is done by a combination of 2/3-rule truncation 
and phase-shifting [16]. As a criterion for sufficient spec- 
tral resolution k max i] > 1 is used [1] with Kolmogorov 
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length scale r\ = (^ 3 /e) 1 / 4 and k max = V2N/3. The as- 
pect ratio L x : L y : L z = 2ir : ir : 2ir was resolved with 
N x (N/2 + l)xN grid points where N = 256 and 512, 
respectively . 

We take periodic boundary conditions in streamwisc 
and spanwise directions, respectively. In the shear direc- 
tion, free-slip boundary conditions are applied, 



du x du z 
u y = 0, and — — = — — = 0. 

ay oy 



(8) 



The free-slip boundary conditions will cause slight devi- 
ations from the transverse homogeneity, an effect which 
decreases with growing Reynolds number [20]. In con- 
trast to the no-slip case, no energy is transfered into the 
flow via the free-slip boundaries and thus an additional 
volume forcing F has to be applied (cf. eq. (3)) in order 
to sustain turbulence. 



space might be nonanalytic in the real space and will vary 
in space and time. 

Second, the grid forcing at small wavenumbers is ap- 
plied that injects a certain amount of energy per time 
unit into the flow given by the rate ej„, [21] 



F#(t) 



u k (t) 



"£ k » e *>k«(*)l s 

Consequently, it would follow 



<5k.k" • 



(12) 



(13) 



for the statistically stationary balance of the turbulent 
kinetic energy when S = 0. The wavenumbers of the set 
K do not coincide with the ones for the shear forcing, 
but are also small. 



Forcing scheme 

A compact formulation of the Navier-Stokes equations 
in Fourier space for our problem is given by 

= N ik (u jv t;Re) + F ik (t) , (9) 

where contains all nonlinear mode coupling contri- 
butions and the dissipative part, — -Re _1 fc 2 Uik(i)- In the 
following we will combine two kinds of forcing, a volume 
forcing that sustains the linear mean flow profile which 
will be referred to as shear forcing and a homogeneous 
forcing that mimics the injection of additional fluctua- 
tions by the grid which will be referred to as grid forcing. 

A handful of Fourier modes is driven now in order to 
give an almost linear profile for (u x ) a,t, he. the forcing is 
with respect to the streamwise velocity component only. 
The driven modes are for the wave vectors k' = (0, 2n + 
1,0) for n = to 5. They form the almost linear mean 
profile, 

, . , . 2a/ 7r\ 8a v-^ cos[(2n + l)y] , . 

= v (y - 2) M - ^ £ L + if ^ 

n=Q v ; 

for y E [0, 7r]. This results in a mean shear rate S = 
2a /it. The shear forcing is chosen in such a way that just 
these modes (which have real parts only for symmetry 
reasons) are held fixed to the amplitude values as given 
by Eq. (10), i.e. 

= ° F^,(t) = -N xU ,(u M ,t;Re). (11) 

These coefficients are set at the begining and do not vi- 
olate the boundary conditions and the divergence-zero 
condition because the driven modes u^k' (t) are functions 
of k' only. The particular kind of forcing in the Fourier 



STATISTICAL STATIONARITY 
Difference to the remeshing method 

The Reynolds decomposition of eqns. (3) and (4) for 
all fields results in 

(u^iuj = -di(p) + ^<«i> + {Fi) - djiviVjm 

where we considered already the time independence of 
the mean profile. At this point, a difference to the case 
with remeshing becomes obvious. All terms of the equa- 
tion will be zero independently of each other due to ex- 
act homogeneity for the latter case. The violation of the 
transverse homogeneity will cause relics in the Reynolds 
balance for the present model and results in 

0=(F^) AiT -dy(v x v v ) A ,T. (15) 

Figure 2 illustrates this behavior. We show the verti- 
cal profile of the shear stress, (v x v v )a,t (upper panel) 
and verify that Eq. (15) is satisfied to a good approxima- 
tion except very close to the boundaries (lower panel). 
The upper panel shows also that our system becomes 
nearly homogeneous only with respect to the transverse 
direction. To conclude, there is a difference between both 
methods for the large scale balance. The question is now: 
how are small-scale statistical properties of the fluctu- 
ating quantities affected in the present approach? As 
demonstrated in Refs. 10 and 11 small scale statistical 
properties such as Reynolds stress magnitudes or higher 
order normalized moments of the transverse derivative 
and spanwise vorticity, respectively, agree with previous 
homogeneous shear flow simulations [9, 17] and experi- 
ments for the lowest Reynolds numbers [5]. 

The forcing keeps the flow in a statistically station- 
ary state with moderate fluctuations of the kinetic en- 
ergy and enstrophy, respectively (see Ref. 11 for a more 
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FIG. 2: Upper panel: Reynolds shear stress — {v x v v )a,t as a 
function of y. Average was taken with respect to time and 
over x — z planes at fixed y. Lower panel: Verification of 

relation (f5). It can be seen that (Fi s ')A,T and d y (v x v y ) a,t 
collapse for almost all y. 
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FIG. 3: Statistical stationarity of the shear flow. Left column 
is for run IV, the right column is for run lie (see Tab.l). Up- 
per panels: Temporal evolution of turbulent kinetic energy 
Ekin(t), energy dissipation rate e(t) (thin solid line), and pro- 
duction of turbulent kinetic energy P(t) (thick solid line). 
Mid panels: Temporal evolution of the ratio P(t)/e(t). The 
thin solid line is the temporal mean of the data and very close 
to one. Lower panels show the temporal variation of the shear 
parameter S*. Temporal means are again indicated by thin 
solid lines. 
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FIG. 4: Temporal evolution of the dimensionless shear pa- 
rameter in the initial phase of the simulation for run IV. The 
run started with isotropic turbulence decaying with a fc~ 5 / 3 
spectral law and was asserted to the constant shear with rate 
S. For comparison we added initial data from Rogers et al. 
[26] (run C128V). 

detailed investigation). In order to demonstrate this 
property of the current numerical scheme, we show 
time series of the turbulent kinetic energy, E k i n (t) — 
q 2 {t)/2 = (v 2 )y/2, the production of turbulent kinetic 
energy, P{t) = —(v x v y )vS, and the energy dissipation 
rate, e(t) = v{(diVj) 2 )y for two runs in Fig. 3 (run lie 
and IV of Tab. 1). The fluctuations of the kinetic en- 
ergy are by about a factor of 2 smaller than those re- 
ported in the long-time remeshing runs. [9, 12] With ge — 
\J((E km {t) - {E kln ) T ) 2 ) T we get a E /(E kln ) T = 24% for 
run IV and 33% for run lie. This ratio decreases even 
further for runs la to Ie. The resulting global energy 
balance for the fluctuating part is then 

fty « -e(t) + P(t) + tin ■ (16) 

It was verified that contributions from shear driving are 
subdominant as well as the transverse flux contributions 
from third order terms arising in the balance. The re- 
duced fluctuations of the kinetic energy might be caused 
due to the free-slip boundaries. As mentioned above, 
small boundary layers do form that cause slight trans- 
verse fluxes into the bulk. [22, 23] 

As can be seen in Fig. 3, the ratio of P/e varies in time, 
but its temporal mean is almost exactly unity and 

(P(t)) T = (c(i)) T - e m , (17) 

holds. We note that this property differs from the 
remeshing simulations where the ratio P/e is constant 
but larger unity. Quantities such as the turbulent kinetic 
energy will grow then exponentially with time which is 
a consequence of (16), i.e. q 2 (t) <~ cxp[(P/e — l)t\. One 
might conclude that the violations of the transverse ho- 
mogeneity cause this ratio to vary around unit value and 
thus to assure statistical stationarity. This point should 
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not be mixed with the exclusion of statistical station- 
arity and downstream homogeneity which is relevant in 
experiments [24, 25] but not in DNS. 

It is interesting to take a closer look at the initial phase 
of the evolution. In Fig. 4 we plotted the initial evolution 
of the shear parameter S* and for comparison we added 
DNS data by Rogers et al. [26]. Both simulations show 
that for the initial phase an (almost) linear growth can be 
detected. It clearly illustrates the non-normal amplifica- 
tion mechanism acting on the isotropic turbulent "back- 
ground" . The streamwise velocity fluctuations grow due 
to the lift-up of streamwise streaks, but will leave the 
energy dissipation rate which is dominant at small scales 
close to its initial magnitude for a while. In Fig. 5 we 
show the streamwise turbulent velocity at two instants, 
an initial snapshot (upper panel) where the fluctuations 
are still isotropic and a later snapshot (lower panel) where 
streamwise streaks have formed. A more detailed analy- 
sis of the regeneration cycle of these streamwise streaks 
and vortices [27] in a shear flow with free-slip boundaries 
was discussed in Ref. 28 and was found to agree qualita- 
tively with other shear flows. 

Integral scale 

In Fig. 6 we show the integral length scale which is 
defined as [1] 

with the energy spectrum, E(k,i). In comparison to the 
simulations that use the remeshing method of Rogallo, it 
can be seen that Ln(t) remains well below the box width 
for our method and does not grow exponentially in time. 
It could be shown recently that strong fluctuations of the 
kinetic energy and the enstrophy, respectively, are related 
to Ln(i) reaching the size of the box length [18]. The 
present studies give L u (t) at about 33% of L y = n and 
are smaller as the reported 72% to 80% of L = 2n.[9] 

INVESTIGATIONS IN THE TWO-PARAMETER 
PLANE 

Estimates for the behavior of S* versus Re 

We turn now to investigations in the parameter plane 
that is spanned by the two essential dimcnsionless param- 
eters in a homogeneous shear flow, the shear parameter 
S* and the Reynolds number Re. Every point (S*,Re) 
in this plane will be denoted as an operating point of the 
particular simulation or measurement, it will drift across 
the plane for the non-stationary case. As was demon- 
strated in Ref. 14, the collected data were rather scat- 
tered over such a plane. Two limiting regimes for the ho- 




FIG. 5: Isosurface plots of the streamwise turbulent velocity 
for the initial phase of run IV. Level sets with opposite sign are 
shown in each picture (grey is positive and black is negative). 
Upper panel: St=0.4. Lower panel: St=9.6. 

mogeneous shear flow could be identified. One is the large 
shear limit which will approach the (linear) rapid distor- 
tion case for S — > oo. A condition for the large shear case 
follows from starting with \vid(uj)/dxi\ » \vidvj/dxi\, 
and it is given by 

S*Re- 1 / 2 > 1 , (19) 

in terms of the large scale Reynolds number Re. The 
second case is the local isotropy limit with the condition of 
sufficient separation of large and small time scales St v <C 
1 . t v = \Jvje is the Kolmogorov time and it follows 

S*Re- 1/2 < 1 . (20) 

The idea is to identify a more systematic dependence 
S*(Re) in the plane, similar to a phase diagram in ther- 
modynamics. It is clear from the beginning that the 
range of parameters that we can cover with the DNS ex- 
periments will be limited, but at least in a certain range 
we can conduct systematic studies. 
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FIG. 6: Temporal evolution of the integral length scale in (t). 
For comparison, the width of the simulation domain in shear 
direction y is indicated by the vertical double-headed arrow. 
The temporal average of in (£) is plotted as a solid line. Data 
are for run IV (upper panel) and run la (lower panel). 



Before doing so, we give simple arguments of which 
functional dependencies (1) and (2) might appear. Dis- 
sipation of energy is due to shear effects and due to grid 
forcing and thus we can set the following relation 



dv- 



Ll 



Li, 



(21) 



where C\ and C 2 are dimensionless constants. The two 
terms on the r.h.s. of (21) model the crossover of the 
energy dissipation rate from a weakly turbulent state to 
fully developed turbulence which is present for moder- 
ate Reynolds numbers. [29] While vv^ ms /Ly determines 
the dissipation at smaller Reynolds numbers, the term 
v^ ms /L y takes over for larger Reynolds numbers (for 
more details see also Ref. 30). The root mean square 
velocity will be determined by both large scale driving 
mechanisms and thus 



= C 3 (L 



,1/3 



(22) 



can be taken by dimensional arguments with yet un- 
known dimensionless constants C3 and C4, respectively. 
Equations (21) and (22) are inserted into the relations 
for the Reynolds number (5) and the shear parameter 
(7), respectively. For both parameters follow functions 
of S and e in in the presence of constants v and L y (cf. 
Eqns. (1) and (2)) 



S* = 



Re = 



LIS 



Gxv + C 2 C 3 L y (L y e m + C A L y S 3 ) 1/3 ' 

C 3 L y (Ly€i n + C^LyS^ 



(23) 



(24) 



One point becomes obvious immediately: there seems 
to be a variety of pathways to explore the S*-Re plane 
and the operating points can be scattered over a large 
fraction. In fact, that is what the collected data points 
showed. [14] 

For simplicity, we study the behavior along curves 
5'=const and ei„=const. The Reynolds number grows 
to infinity for both cases because one of the sources for 
an increase of turbulent fluctuations is always present. It 
reduces to 



Re = 



r 7-4/3 1/3 



(25) 



for turbulence without shear and marks the origin of 
curves S*(Re) in the plane (cf. relations (19) and (20)). 
We get the following limits for the shear parameter, 



lim S* 

S=const, €i n — »oo 



0, 



lim 

€in = const, S— > 00 



s* 



(26) 
(27) 



In limit (26) we observe that the shear parameter goes to 
zero with a power C(e irl 1/3 ). Physically, this means that 
the grid driving becomes more and more important while 
shear effects decrease, i.e. the system approaches the 
isotropic limit S*(Re — i?e )~ 1//2 <C 1, which is the case 
to a good approximation in grid turbulence. [14] On the 
other hand, when keeping ei n fixed but increasing exter- 
nal shear rate S, the estimate yields a non-zero constant 
value if there is no further Reynolds number dependence 
hidden in the four prefactors Cj. This means that for a 
statistically stationary and nearly homogeneous system 
an asymptotic value of S* is found. 



Simulation results 

The outlined estimates are tested by numerical ex- 
periments. In Table I, the different parameter sets and 
the quantities are summarized that are necessary for the 
calculation of S* and Re. The DNS runs are ordered 
in three different series, No. I varies along the isoline 
(S = const, ei n ), while series II and III run along the 
iso-parameter line (S,e in = const.). Such investigation 
becomes rather expensive because every single data point 
in the parameter plane is one long-time DNS run for a 
statistically stationary shear flow. 

In Fig. 7, we have collected the operating points 
(S* ,Re) of the runs and added the error bars. The inset 
illustrates the corresponding variations in terms of the 
outer parameters. We fitted the dimensional estimates 
to the whole set of data points. The dashed and the 
dotted lines are the fit results of (23) and (24) to series 
I, II, and III. The additional lines illustrate how points 
(S, ei n ) are mapped to (S* , Re) with the fitted constants. 
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Run No. 


s 


tin 


Vrms 


e 


Re 


S* 




v N 


10r 


Ia 




1/tt 


0.010 


0.47 


0.017 


444 


4.2 


4.64 


1/300 256 




lb 




1/tt 


0.025 


0.55 


0.032 


520 


3.0 


3.96 


1/300 256 


9 


Ic 




1/tt 


0.040 


0.66 


0.051 


622 


2.7 


3.52 


1/300 256 


8 


Id 




1/tt 


0.055 


0.67 


0.063 


629 


2.2 


3.34 


1/300 256 




Ie 




1/tt 


0.100 


0.86 


0.114 


810 


2.1 


2.88 


1/300 256 


7- 


Ila 







0.010 


0.34 


0.010 


325 


0.0 


5.29 


1/300 256 


6 


lib ( 


=Ia) 


1/tt 


0.010 


0.47 


0.017 


444 


4.2 


4.64 


1/300 256 




lie 




3/tt 


0.010 


0.94 


0.111 


891 


7.7 


2.90 


1/300 256 


*co 5 - 


lid 




6/tt 


0.010 


1.49 


0.480 


1404 


8.8 


2.01 


1/300 256 


A - 


He 




8/tt 


0.010 


3.39 


2.998 


3197 


9.8 


2.54 


1/300 512 




Ilia 







0.025 


0.48 


0.025 


457 


0.0 


4.21 


1/300 256 


3 


Illb 


;=ib) 1/tt 


0.025 


0.55 


0.032 


520 


3.0 


3.96 


1/300 256 




IIIc 




6/tt 


0.025 


1.31 


0.393 


1235 


8.4 


2.11 


1/300 256 


2 


IV 




1/tt 





0.36 


0.005 


452 


8.2 


5.07 


1/400 256 


1 - 



TABLE I: Parameters of the numerical experiments. For con- 
vinience we listed the runs in three series indicating different 
pathways in the parameter plane. Run IV was an additional 
long-time run for investigations on statistical stationarity. 



It can be seen that the lines display the trends of the 
data points quite well, but do not match perfectly with 
the data points. Clearly, the data base is rather sparse 
such that for a multi-dimensional fit problems can arise. 
Second, we do observe that the error bars grow in mag- 
nitude for the larger shear rates when the fluctuations 
of large scale velocity tend to grow because the effect of 
non-normal streak lift-up gets more pronounced. 

The numerical experiments support a saturation of 
the shear parameter for larger Reynolds numbers which 
would mean that a flow with a given injection rate can 
never cross the large shear limit boundary of S*(Re — 
Reo)^ 1 ^ 2 S> 1 when S is increased. Clearly, at the present 
stage this point cannot be fully resolved, but we do ob- 
serve a saturation along II and III. It might be the case 
that even for the present system the finite size of the 
simulation domain affects the results eventually although 
the integral scale remains well below the box size. On the 
other hand and to our knowledge, all DNS operated in a 
range of the shear parameter of S* < 10 that could be 
exceeded only transiently when the system was strongly 
non-stationary and therefore close to the rapid distortion 
case [19, 26]. 

Further support for the ansatz (22) is given by the 
following fact. When inserting (22) into (21), one gets e ~ 
ej„ + const S 3 to leading order plus further subdominant 
terms in each of the variables. In Fig. 8, we plotted e — ej„ 
over shear rate S for all three series and observe that for 
the larger shear rates a cubic power law fits quite well. 
The scaling is consistent with a direct proportionality 
between S and v rms (cf. (21)). 
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FIG. 7: Operating points of statistically stationary homoge- 
neous shear flow simulations in the Re-S* parameter plane 
when iso-parameter lines were followed as indicated in the 
inset. Asterisks are data points of series I, open circles of 
series II, and open triangles of series III. The error bars 
(±<r/2) are given in the panel. They follow from as* = 
S*(2a Vrm Jv rm3 +er e /e) and a Re = Rea Vrms /v rms , respec- 
tively. The two data points with S* = have error bars 
with respect to Re which are about the size of the symbol. 
The mapping of the isoline mesh of the inset to the S*-Re 
plane is also shown by additional dotted and dashed lines 
lines. The constants are Ci = 1.5, C2 = 0.35, C3 = 1.29, and 
C 4 = 0.015. 




FIG. 8: Double logarithmic plot of e — e in over S for all three 
data sets. The symbols are the same as for Fig. 7. The data 
were fitted (least square) with a function AS 2 + BS 3 where 
A — 10~ 4 and B = 0.15. The dashed line shows e — e in = 
0A5S 3 . 
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FIG. 9: Normalized derivative moments Ms(dv x /dy) (up- 
per row) and M±{dv x / dy) (lower row) as a function of the 
Reynolds number Re (left column) and the shear parameter 
S* (right column). The data from series I and II of Table 1 
are plotted. 



We turn back to the original scope of these investi- 
gations to get a better understanding of the return to 
local isotropy in shear flows. We monitored derivative 
moments along the parameter curves. The normalized 
n-th order moment of the transverse derivative of the 
streamwise turbulent velocity component is given by 



M n {dv x /dy) 



{{dv x /dy) n ) 
{{dv x /dyf)^ 



(28) 



The cases n — 3 and 4 arc studied and the results 
are summarized in Fig. 9. The flatness (n = 4) in 
the lower row grows slowly with Reynolds number and 
shear parameter, respectively. The Reynolds number 
trend is thus in agreement with previous investigations 
[9, 10, 11, 17]. More interesting arc the results for the 
third order moment, the skewness. While it decreases 
strongly with Re along path I it starts to decline slowly 
only for the largest values of Re or S* along path II. We 
note that for both series the Reynolds number is grow- 
ing while their decay differs. A simple projection onto 
Re gives thus different results and underlines the idea 
of taking into account the shear rate as well. This is 
also supported when plotting all skewness data over S*. 
Now, a trend of the data to have larger skewness val- 
ues with growing shear parameter is observed (cf. upper 
right panel of Fig. 9) for the majority of the data set. 

In Ref. 14 the data were discussed with respect to 
the following question: Is a large Reynolds number re- 
ally large if the shear parameter is large as well? Large 
is meant there in the sense of being close to the local 
isotropy limit. Based on our findings and on (27) one 
can draw the conclusion that large means then 



Re » (S^f 



1 



( - / 2 ( - / 3 4 



(29) 



which can be further simplified. For larger Re one would 
expect C*2 ~ 0(1) because this constant is nothing else 
but the dimensionless energy dissipation rate (see ansatz 
(21)). A similar argument should hold for C3. Thus 
Re » I/C4 3 would follow eventually. The remaining 
constant C4 which relates the turbulent fluctuations to 
the shear rate cannot be evaluated by such arguments. 
The present DNS give C4 <C 1 consistently (see the cap- 
tion of Fig. 7). This can be an interesting point to be 
answered when the data base is more comprehensive. 



SUMMARY AND OUTLOOK 

We have presented systematic studies of the relation 
between the shear parameter and the Reynolds number in 
a turbulent nearly homogeneous shear flow. The present 
system allows for investigations of a statistically station- 
ary flow with mild fluctuations of the turbulent kinetic 
energy. The ratio of the turbulence production to the en- 
ergy dissipation was found to vary around unity and the 
integral length scale remains well below the box exten- 
sion. Two parameters, the shear rate S and the energy 
injection rate ej n , were varied by keeping one constant in 
each case. We monitored the resulting shear parameter 
S* and Reynolds number Re. It was found that they 
evolve on different pathways across the two-parameter 
plane. Dimensional estimates are in qualitative agree- 
ment with the results of the DNS. They suggest that an 
asymptotic shear parameter value is reached for the sta- 
tistically stationary case. The reason for such saturation 
might be due to the finiteness of the considered volume, 
namely that the linear mean profile in a homogeneous 
shear flow causes always an outer flow gradient scale that 
is beyond the box extensions (strictly speaking this scale 
is infinite). For any finite system the initially observed 
self-similar growth of quantities with time is interrupted 
once Ln comes closer to the box size and consequently 
the shear parameter settles to a finite value. 

Another way of such a two-parameter discussion is to 
use the shear parameter S* and the Corrsin parameter 
S c = St v [31] where the Taylor microscale Reynolds num- 
ber follows to R\ ~ S*/S c . [12] Clearly, our results can 
be transformed into the latter frame. The large shear 
limit from Ref. 14, S*R^ 3> 1, would thus translate to 
S c 3> 1 and the local isotropy limit, S*R^ K 1 <C 1, follows 
to 5 C <C 1. The latter inequality states then that the 
shear time scale S^ 1 and the dissipative time scale t v 
have to be separated far enough of each other. 

It was shown that the third order derivative moments 
decay with different trends with respect to the Reynolds 
number when they are monitored along different path- 
ways in the two-parameter plane. There is practically 
no decay over a wide range when the same data are 
shown with respect to the the shear parameter. All this 
might explain the different interpretations of measured 
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moments in terms of a return to local isotropy. 

Clearly, the present the Reynolds numbers are mostly 
small and the data base is rather sparse. It has to become 
more comprehensive for some reasons: first, one would 
like to see of these trends persist to higher Reynolds num- 
bers, especially the behavior of the shear parameter S*. 
Secondly, the four parameters C\ to C\ are assumed to 
be constant. We do not know if a Reynolds number de- 
pendence is hidden there; the current data did not allow 
for drawing conclusions on that. Thus, we consider this 
investigation as a starting point. Further extensions be- 
side these two points are possible. The so-called grid 
forcing scheme gives us a tool into hands which can be 
extended to cases where spatial and temporal patterns 
are used for the excitation of turbulent fluctuations, an 
approach that recently started in terms of active grids in 
wind tunnels. 
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